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We study the electronic transport through a spin-1 molecule in which mechanical stretching 
produces a magnetic anisotropy. In this type of device, a vibron mode along the stretching axis 
will couple naturally to the molecular spin. We consider a single molecular vibrational mode and 
find that the electron-vibron interaction induces an effective correction to the magnetic anisotropy 
that shifts the ground state of the device toward a non-Fermi liquid phase. A transition into a 
Fermi liquid phase could then be achieved, by means of mechanical stretching, passing through an 
underscreened spin-1 Kondo regime. We present numerical renormalization group results for the 
differential conductance, the spectral density, and the magnetic susceptibility across the transition. 



I. INTRODUCTION 

The electronic properties of nanostructures depend 
critically on their symmetries. As a consequence, the 
ability to modify these microscopic symmetries makes 
it possible to drive the system through different physi- 
cal regimes at will, possibly resulting in quantum phase 
transitions (QPTs) as the system visits different ground 
states.El 

Quantum dot^^^ and molecular devices consist- 
ing of cornplex molecules deposited on metallic break 
junctionsSHSl are good examples of systems where phys- 
ical regimes can be explored by tuning their parame- 
ters. The electronic transport through molecular devices 
is very sensitive to the hybridization of the molecular en- 
ergy levels with the bands of the metallic leads in the 
break junction, as well as to electron-electron (e-e) and 
electron-vibron interactions within the molecule.'^i'lSI 

Success has been achieved in the past few years in 
controlling the properties of molecular junctions through 
tuning of the molecular levels by different means. It has 
been shown, for example, that a gate voltage is capable of 
inducing a singlet-triplet transition in a Cgo (buckyball) 
molecule trapped between metallic leads, due to distinct 
hybridizations of these states with the electronic states 
of the leads.l^ 

In recent experiments by Parks et al.^ it was 
shown that mechanical stretching of the spin-1 molecule 
Co(tpy-SH)2 (4'-mercapto-2, 2':6',2"-terpyridine) along 
the transport axis in a break junction setup, can be used 
to control the magnetic properties of the molecule. The 
stretching induces a splitting of the spin-f triplet ground 
state,^^"^!] raising the energy of the doublet with spin pro- 
jection Sz = ±f by as much as 4 meV with respect to the 
state of Sz = and leaving the latter as the molecular 
ground state. 

This magnetic anisotropy is critical to the low- 
temperature transport through the molecular junction in 
this experiment. In the absence of anisotropy, the sys- 



tem exhibits an underscreened spin-1 Kondo effect, sig- 
naled by enhanced conductivity for temperatures below 
the Kondo temperature T^.[22l Any positive (hard-axis) 
anisotropy breaks the ground state degeneracy of the iso- 
lated molecule and drives the system into a Ferm i liqu id 
ground state with an associated low conductance.'^SSSl 

The converse case of negative (easy-axis) anisotropy 
could be reached, in principle, by compressing the 
molecule. This situation would set the S'^ = ±f doublet 
as the molecular ground state, yielding non-Fermi-liquid 
(NFL) behavior.24 26-28 

In this context, visiting different ground states be- 
comes a matter of adjusting the magnetic anisotropy. 
The system will go from an ordinary Fermi liquid (FL) 
in the hard-axis regime to a NFL in the easy-axis regime, 
passing through an underscreened Kondo ground stated 
— a singular Fermi liquid (SFL) — when the triplet is 
exactly degenerate.l^ These different regimes and the as- 
sociated QPTs can be studied as function of a single pa- 
rameter: the magnetic anisotropy. 

All of the effects described above can be caused by 
static deformations of the molecule. Moreover, it is to be 
expected that dynamical effects may arise via an analo- 
gous coupling between the spin and the mechanical de- 
grees of freedom of the moleculcI^SI As we show below, a 
coupling between molecular vibrations and spin will in- 
duce a deformation in the molecular ground state, open- 
ing access to an easy-axis regime. With this motivation, 
in this paper we study a model that encompasses the 
anisotropy regimes described above, and considers in ad- 
dition the mechanical degrees of freedom of the molecule 
through a vibrational mode. In the same way as the 
static deformation, vibrations along the axis couple nat- 
urally to the spin projection of the molecule. In the case 
of an isolated molecule, we find that the coupling to the 
vibrational mode indeed opens access to the easy-axis 
regime, i.e., the spin-f doublet of S'2 = ±f becomes the 
ground state. When the molecule is coupled to leads, we 
find that the magnetic anisotropy is further renormal- 
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FIG. 1. (Color online) Phase diagram of a deformable spin- 
1 molecular break junction in the presence of a local vibra- 
tional mode. A rich variety of phases can be found at tem- 
peratures below Tk, at which Kondo physics dominates the 
transport properties of the system. The effective anisotropy 
Aefs, induced by both static and dynamic deformations of 
the molecule, drives the system through different quantum 
phases. 



ized. The resulting effective anisotropy could be tuned 
to explore a variety of ground states which we will dis- 
cuss below. We studied the different highly-correlated 
regimes of the model by means of numerical renormal- 
ization group (NRG) calculations .^^^ll^ll The results are 
summarized in the phase diagram shown in Fig. [T] to 
which we will come back later. 



II. MODEL 

We model the molecular device shown in Fig. [2ja) as 
depicted in[2jb). The Hamiltonian of the spin-1 molecule 
can be represented by a 2-orbital model of the form 



Ho 



{etii + Uni-^Hi^) ~ J Sa- Sb- 



(1) 



Here, i = a,h are the two degenerate molecular orbitals 
with energy e and intra-orbital Coulomb repulsion [/, 
jT-io- = dl^dia- and dl^{dicr) the creation (annihilation) op- 
erator of the corresponding orbital; rii = ni^+rii^, and Si 
is the spin operator asociated to orbital i. For simplicity, 
we consider the electron-hole (e-h) symmetric case where 
£ = —U/2 and the Fermi level of the leads is £f = 0. The 
ferromagnetic (J > 0) coupling enforces Hund's rule, set- 
ting the spin-1 triplet as the ground state of the molecule. 
These states are defined as 



|T, +l) = ltatb> 



\T, -1) = \iaU). 



\T,0) = ^(^\taU) + \iaW), 

where \(Ta (Jb) are the states with one electron on each 
orbital. 




\A[z) St 



= 0, ± 1 



- 

-1,+1= 



-l.+l 





^eff 



FIG. 2. (Color online) Spin-1, deformable molecular device, 
(a) A break junction supports a spin-1 molecule that can be 
stretched by mechanical means, (b) The molecule is modeled 
by two coupled orbitals, only one of them connected to metal- 
lic leads. The spin-1 regime is enforced by an inter-orbital fer- 
romagnetic coupling. Magnetic anisotropy, induced by static 
and dynamic (vibrational) stretching of the molecule, is ac- 
counted for by the term A{z)Sl in the Hamiltonian. (c) Low- 
est lying level energies of the isolated molecule, obtained for 
different values of the effective static anisotropy. 



The deformation-induced anisotropy arising from the 
breaking of the octahedral symmetry can be written as a 
function of the elongation along the z axis, as A{z)Sl. A 
purely static model [A{z) = Aq] of this kind was applied 
in reference [3T] to describe the system of reference [S]. 
For small oscillations around the equilibrium position zq 
we take A{z) = [Aq + Ax 5z) Si, where 5z = z — zq. The 
coupling Aq represents the static deformation, whereas 
Ai arises due to the e-ph coupling. We choose units such 
that 5z = a-\- a\ where a and are phonon operators. 
The Hamiltonian of the molecule including the electron- 
vibron interaction becomes 

Hm = Ho + Ao Sl + Ai Sl (a -f a)) + a) a, (3) 

with LJa the phonon frequency. By means of the unitary 
transformation 



Z^ = cxp|-^5,2(a-at)|, (4) 
we diagonalize (|3|, which adopts the form 



Hm = Ho + AS'l + loq 



(5) 



where b = a + {Ai/ijJq) S^ is a displaced phonon opera- 
tor. Note that the magnetic anisotropy of the isolated 
molecule 



A = Aa- Al/ujo, 



(6) 



is negative in the absence of a static distortion {Aq = 0) . 
The lowest energy eigenstates of ^ for each molecular 
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spin projection are 

|T, ±1; 0) = IT, ±1) |0) , |T, 0; 0) = |T, 0) |0) , (7) 
where \n) is an eigenvector of the operator a^a, and \h) — 



e"o 



i(a-an 

n ^ ' 



i) is an eigenvector of b^b, with eigenvalue 



We now leave the isolated molecule (molecular limit) to 
explore the consequences of coupling the metallic leads to 
the molecular orbital a, through the hybridization term 



(8) 



This "hanging-level" arrangement correctly describes the 
low-energy behavior of the system, and other configura- 
tions can be related to it by means of a level rotation, so 
there is no loss of generalityPil Because we assume iden- 
tical right (R) and left (L) leads and couplings, we have 
defined the lead-symmetric operators 



-Rfe, 



(9) 



which are the only combinations that couple to the 
molecule. Their anti-symmetric counterparts contribute 
only a constant energy term to the Hamiltonian. 

In the displaced basis the full transformed Hamiltonian 
reads 



H = Hn + ASt 



- Wo 



b% 



E 



ka 



H.c. 



(10) 



with a hybridization operator 



/2F exp ( ^ (a - a)) [l - 2nas + 4a5,^] 



Although the expression for the transformed Hamilton- 
ain is more complicated than the starting one, we can 
appreciate an interesting feature. Namely, that it cou- 
ples the ±1 and spin projections of the triplet to the 
electronic states of the band with different strengths. As 
we will see below, this asymmetry leads to a contribution 
to the total magnetic anisotropy. 



A. Analytical results 

To expose the splitting of the triplet due to the spin de- 
pendent couplings to the band, we perform a Schrieffer- 
Wolff transformatiorPSl in the limit of large U and J. 
Here we set A — Q so that the isolated molecule ground 
state is triply degenerate. We obtain an anisotropic 
Kondo HamiltoniarP^ 

Hk = A^Sl + 4s, + (s, + sy Sy) , (12) 

with exchange couplings 



J 



K 



E 



4V^\(b\ n)\ , n , ^ 

IV I n T-L _ „ (.0 j /V^, (13) 



u 

2 



where = AV'^/{^ -I- ^) is the Kondo coupling in the 
absence of e-ph coupling {Ai = 0), and 



1 V I (01 



(0| n)| 



u 

2 ' 4 



riLOo 



> 0, (14) 



where 



|(6| n)\' = 



Wo 



(15) 



The anisotropy term Ad is a consequence of the hy- 
bridization shifting the state \T, 0; 0) down in energy 
further than the states \T, ±1; 0), for any state of the 
band. 

The e-ph interaction induces a reduction of both 
Kondo couplings that become anisotropic {J^^Jk < 
J^). In the experimentaly relevant case C/ 3> we 

have < jI^. For weak e-ph interaction (Ai/wq <C 1) 
and U ^ we obtain: 



Ai 



1 - 



Wo 



Wq 

M 

Wo 



(16) 
(17) 
(18) 



The expresions for and Ad are also valid in the strong 
e-ph coupling regime {Ai 3> wq), provided A^/wq <C U, 
while is exponentially supressed [see Eq. ( 13 )] making 
the Kondo couplings strongly anisotropic. 

As a result of the correction A^, the degeneracy of 
the triplet is broken preventing the underscreened spin-1 
Kondo effect from taking place at A = 0. Therefore, a 
shift in the static anisotropy will be needed to restore 
the full spin-1 degeneracy. It is also important to point 
out that the anisotropy of the Kondo couplings leads to 
effects similar to those obtained by a static magnetic 
anisotropy term (Ao). In our case of > > 0, 
this a nisotr opy results in a shift toward the easy-axis 
regime.™! 

To summarize, the exact solution of the isolated spin- 
1 molecule, along with a perturbative analysis to second 
order in the hybridization, suggest that the main effect of 
the electron-phonon interaction considered in this break 
junction setup is to reduce the Kondo couplings, and to 
introduce corrections to the magnetic anisotropy. The 
different contributions to the magnetic anisotropy result 
in an effective anisotropy given by 

Aoff = Ao - A^/wo + Ad + (5A, 

where 5 A stems from the anisotropy of the Kondo cou- 
plings. 
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2.5, with no fewer than 1200 states 
and a constant hybridization F = 7r|y|^/2Z? — 



III. NUMERICAL RESULTS 

We analyze the different magnetic anisotropy regimes 
of the system as a function of Aq, Ax and Wq, by means 
of NRG calculations. In this non-perturbative method, 
the continuum of electronic states in the leads is logarith- 
mically discretized according to a parameter A > 1. It 
is then mapped onto a chain of states with exponentially 
decaying hopping terms that is diagonalized iteratively, 
defining a renormalization group transformation to the 
low energy spectrum. '^^ 

All calculations shown below were obtained with dis- 
cretization factor A 
kept 

0.005 1?, where D is the half-bandwidth. We choose the 
parameters e = —11/2 — — 5r, with ferromagnetic cou- 
pling set to J = 0.2 r. Our conclusions, however, remain 
qualitatively unchanged for other parameters within the 
Kondo regime (such as J > F). The spectral densities 
were calculated following Ref. [57] . 

For Ai^ = Ai = 0, the system is expected to have 
an underscreened spin-l Kondo ground state; we have 
verified this with NRG calculations. This regime serves 
as a reference for comparison with our results in other 
regions of parameter space. 

The infinite-dimensional boson sector of the Hilbert 
space has to be truncated for the numerical calculation. 
We analyze the regime of wq ^ and Ax/uq < 1, so 
that relevant phonon excitations include transitions from 
the ground state with boson component given by |0) to 
states with occupation n; these have amplitudes given 
by Eq. (15). This Poisson distribution peaks at n = 




0.30 
0.20 
<^ 0.10 
0.0 



Ai/ujq < 1 and then falls to zero rapidly, suggesting a 
cutoff at phonon occupation n ~ 1.^^ 



FIG. 3. (Color online) (a) Spectral density at orbital a for 
Aq = Ai = (reference system). The screening of the 
molecular triplet ground state by the band is characterized 
by the sharp Kondo resonance of width at the Fermi level 
{uj — 0), followed by shoulders at approximately the singlet- 
triplet (red dashed lines), and the charge (blue dashed lines 
in the inset) excitation energies of the isolated molecule. The 
Kondo effect yields unitary conductance at low (^ T^) tem- 
peratures, (b) Same as (a), but for finite magnetic anisotropy 
(^0 = tjo = 5Tff, Ai from 2T^ to 11 ). In the easy 
axis regime of AaB < (red curves), the system presents a fi- 
nite amplitude of the spectral density at the Fermi level. The 
hard-axis regime {AcS > 0; blue curves) presents a dip at the 
Fermi level, keeping the system in a Coulomb blockade. The 
underscreened Kondo state is recovered when AcS — (black 
curve). 



A. Transport properties at low temperatures 

The low-temperature, linear transport properties of 
the system can be obtained in terms of the spectral den- 
sity at orbital a, Pa( ui) = Pai^)^ which relates to the 
conductance througlP^lifll 



G{T) 
Go 



ttF 



duj 



(19) 



where Gq = 2e^ /h, and /(w, T) is the Fermi distribution. 
In the reference system, where Aq = Ai — 0, Pa(w) shows 
a sharp Kondo resonance of width T]^ ^ 10~^ D at the 
Fermi level^'^'^ (e^^-, set here to zero) as shown in Fig. 
|3]ja). At low temperatures the conductance reaches the 
unitary limit. Shoulders in the reference spectral density 
appear at approximately the singlet-triplet excitation en- 
ergy, UJ ^ ± J, and peaks at the charge excitation energies 
w - ±U/2. 

With this in mind we move on to analyze Fig. [sjjb) in 
the hard-axis regime {A^g > 0). The spectral density has 
two peaks close to the Fermi level separated by a sharp 



dip, which according to ( 19 ) puts the device in a trans- 



port blockade with a vanishing zero-bias conductance at 

In the easy-axis regime {Acs < 0), conduction elec- 
trons will scatter off the molecule's ground state,'23 the 
doublet |T, ±1). At zero temperature Paisp) presents a 
non-zero amplitude, which means that the system will 
conduct at zero bias, although not unitarily as in the 
case of zero anisotropy. The different regimes described 
above are shown in Fig. |3jb) from NRG calculations in 
which Aos was tuned by varying Ai, at constant static 
anisotropy Aq and phonon frequency ojq. 

Curves of conductance as a function of temperature, 
corresponding to the parameters of Fig. [sjb), are shown 
in Fig. |4]ja), whereas in Fig. |4]jb) the same regimes are 
explored by varying Aq instead. As mentioned above, 
tuning Aq in molecular devices has been achieved exper- 
imentally by Parks et aL,'^ through stretching. Assuming 
that the value of Aq increases from zero by stretching the 
molecule, it could be used to change the sign of the effec- 
tive magnetic anisotropy. This is depicted in Fig. Qb) , 
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FIG. 4. (Color online) Conductance as a function of tempera- 
ture, for different values of AcS- In (a), varying A\ from 2T% 
to 12 for ujQ = Aq = ^T%, and in (b) by varying Aq from 
4,T% to 9ri for = cjo = 5r^. In (c), as a function of Aq 
for different isotherms. The zero-temperature conductance 
vanishes for positive anisotropy, due to the lack of available 
states for transport in the molecule at the Fermi level (see 
Fig. [3|. As yleff — >■ the conductance curves vary as indi- 
cated by the blue arrows until the unitary limit is reached, 
signaled in the plot by a plateau that extends to lower tem- 
peratures. Unitarity is then lost due to many-body scatter- 
ing of the leads' electrons off the Sz = ±1 degenerate ground 
state of the molecule at negative anisotropies. As AcS goes 
from zero to negative the conductance curves follow the trend 
indicated by the red arrow. In (c) we show isotherms for 
the conductance as a function of (corresponding to the 
curves in (b)), which can be directly related to experiments 
in which the conductance is measured in a stretchable device, 
at constant temperature. The quantum critical point becomes 
apparent at very low temperatures, where a sudden drop in 
conductance signals the change of sign of AeS- The Kondo 
temperatures of these cases are lower than the reference scale 
by one order of magnitude, as can be seen by the onset of 
the conductance plateaus below T — T%. 



where we begin with a positive value of Aq such that the 
effective anisotropy is negative, and then increase Aq un- 
til we reach the hard-axis regime. Experiments such as 
those of reference [6J may perhaps be more easily related 
to Fig. |4][c), where we show the zero-bias conductance 
at constant temperature, for different values of Aq. The 
sudden drop in conductance is a clear signature of the 
transition from hard-axis to easy-axis behavior at zero 
temperature. 



B. Effective anisotropy and underscreened Kondo 
effect restoration 

As we will see below, the numerical results confirm 
that the main effect of the e-ph coupling considered is to 
renormalize the magnetic anisotropy toward the hard- 
axis regime and reduce the Kondo exchange couplings. 
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FIG. 5. (Color online) Effective magnetic-moment-squared, 
/i^ = T x{T), as a function of temperature for different values 
of the magnetic anisotropy. Solid lines are for a system with 
e-ph interaction (^i =^ 0), whereas dashed lines correspond to 
a system without e-ph interaction, and with a static magnetic 
anisotropy tuned to match the effective anisotropy of the e- 
ph interacting system. As expected from the e-ph mediated 
reduction of the Kondo exchange couplings, the first-stage 
Kondo temperature is lower in the e-ph interacting system 
[solid arrows in (b)], that is, Tk < T%. This in turn increases 
the second-stage Kondo temperature [arrows in (a)], as 
expected from Eq. (201. 



In fact, the static anisotropy term can be tuned to re- 
cover the underscreened spin-1 Kondo effect. We define 
A = AcQ — A and calculate it numerically through an 
analysis of the effective magnetic-moment -squared, p?, 
at low temperature.-^ The onset of Kondo screening is 
signaled by a drop in the molecule's magnetic moment at 
temperatures below some Tk{Aq, Ai, ujq). Notice that 
Tr- < due to the reduced Kondo couplings [see Eqs. 
(13)]. An example of this is shown in Fig. [sjb). 

As the molecule is coupled to a single conduction elec- 
tron channel, the electrons in the leads are able to screen 
one half of the molecule 's spi n, leaving an asymptoti- 
cally free spin-1/2 objecpSBSI andjaroducing a plateau 



ipr 

^2(rp ^ j^^-j ^ Q pjg [5j Further screening 

associated to a second-stage Kondo effect is known to 
arise when a net positive anisotropy < A^s ^ Tk is 
present. We have H^{T) ~ for temperatures below a 
second stage Kondo temperature T^, given bjPH 



Tk = ciTkG 



(20) 



with constants ci, C2 1 that depend weakly on the 
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parameters. 

On the easy-axis side, the Kondo screening is sup- 
pressed by the energy gap between the molecular ground 
state |T, ±1; 0) and the first excited state |T, 0; 0). In 
other words, Kondo screening becomes increasingly less 
efficient as we go farther away from full degeneracy of the 
triplet, and we obtain 0.25 < fi^ < 1 for all temperatures 
(not shown). 

In order to quantify A^g and estimate A, we carry out 
NRG calculations at fixed Aq, varying Ai and ujq. The 
transition from a hard to an easy axis is observed as a 
sudden jump (at zero temperature) from the hard-axis 
value of /x^ = to /i^ ~ 0.25 in the spin-1 Kondo regime, 
as shown in Fig. [6f^a). 

For the molecular limit where the magnetic anisotropy 
coefficient is exactly given by A, the change of regime is 
indicated by the red curve in Fig.[6][a) and (c), given by 
A = 0. When the coupling to the band is taken into ac- 
count, the transition occurs for Acs = and is signaled 
by ^^(T 0) ~ 0.25. We indicate this transition with 
black dots in the figures, and find A = A^s — A > 
as expected from our discussion above. The values of A 
for every luq are presented in Fig. and (d) for their 
respective maps. Notice that the points at larger ujq cor- 
respond to smaller Ai/ujq, which explains the decreasing 
trend of A with increasing ujq. 
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FIG. 6. (Color online) (a) and (c): Magnetic-moment- 
squared of the molecule at zero temperature (/i^(T — >■ 0)) as a 
function of Ai and loq, for (a) Aq = ST]^ and (c) Ao — 6T^k- 
In each case, the red curve (Ai = y/ujoAo) indicates the pa- 
rameters of full degeneracy of the spin-1 triplet in the molec- 
ular limit. The black dots indicate the parameters where the 
underscreened spin-1 Kondo effect is recovered {fi^ — 0.25). 
(b) and (d): Correction to the magnetic anisotropy due to 
the coupling to the leads for the parameters of (a) and (c), 
respectively. 

The e-ph interactions suppress the Kondo couplings, 
reducing the spin-1 Kondo temperature Tr- of the device. 



This results, according to Eq. ([20j), in an enhancement of 
. This can be observed in Fig. [5j where we compare 
our device to one with only static anisotropy, and find 
that they behave similarly, differing only in the values of 
Tk and T^. 

As further verification, we utilize the values of A ob- 
tained from Fig. [6 to compute values of A^g for sub- 
stitution into Eq. (20), and the degree of agreement is 
excellent, as can be appreciated in Fig. [tJc). 

The excellent fit to the theory reassures us of our 
picture, in which all the effects of the electron-phonon 
(e-ph) interactions can be absorbed into an effective 
anisotropy term, and a reduced hybridization due to po- 
laronic effects. 
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FIG. 7. (Color online) Second-stage Kondo temperature, T^, 
for different values of Ai and uq. (a) A fast fall of as 
a function of Ai/uq is observed, consistent with a vibronic 
induced shift of the magnetic anisotropy and Eq. (201. (b) 
Plotting against Af/ujo — Ao — A, makes the presence 
of a correction A to the anisotropy evident; the dotted lines 
indicate the value of Ai where the molecular limit anisotropy 
A changes sign, for the corresponding values of Aq. The nu- 
merical results are clearly shifted toward positive values of 
^cft. (c) Verification of Eq. (20 1 within our picture of effec- 



tive anisotropy and suppressed Kondo exchange couplings, for 
fixed Ai and too, varying Aq. We substitute A^a in the equa- 
tion, with A computed from numerical results such as those of 
Fig. [6|b) and (d) . The slope of each curve is proportional to 
— with Tk the Kondo temperature corresponding 

to the curve. 



IV. CONCLUSIONS 

We have studied the behavior of a stretchable spin-I 
molecule deposited on a break junction, that presents 
a vibrational mode along the junction axis. We per- 
formed NRG calculations to study the Kondo physics 
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of this system; its signatures in transport and thermo- 
dynamic quantities, as well as the nature of the system's 
ground state. We find that the vibrational degrees of free- 
dom induce a negative magnetic anisotropy that shifts 
the ground state of the undcformed molecule toward an 
easy-axis regime, which would be accessible to the static 
molecule only through compression. Polaronic correc- 
tions arising from the electron-phonon interactions sup- 
press the Kondo exchange couplings of the molecular spin 
with the leads' fermionic states, in a spin-asymmetric 
fashion. This reduces the Kondo temperature of the de- 
vice and introduces an effective correction to the mag- 
netic anisotropy. 

Static stretching of the molecule could then be used to 
explore a QPT from the non-Fermi-liquid ground state of 
the easy-axis regime, to the Fermi liquid of the hard-axis 
regime, visiting an underscreened Kondo critical point at 
zero anisotropy. The resulting phase diagram of the sys- 
tem is presented in Fig. [T] The different phases exhibit 
clear conductance signatures that can be explored exper- 



imentally. 

From our analysis, we may expect the un-stretched de- 
vice in the experiment of Parks et al., to be in an easy- 
axis regime due to the coupling to molecular vibrations. 
In that case, a careful analysis of the low-temperature be- 
havior of the conductance, such as that indicated inf Fig. 
^c), would make it possible to identify the signatures of 
this QPT in a well-controlled environment. 
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